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Abstract 

Motivated by the potential applieations of their intrinsie eross-eoupling prop¬ 
erties, the interest in multiferroie materials has eonstantly inereased reeently, 
leading to signifieant experimental and theoretieal advaneements. From the the¬ 
oretical point of view, recent progresses have allowed to identify different mecha¬ 
nisms responsible for the appearence of ferroelectric polarization coexisting with 
- and coupled to - magnetic properties. This chapter aims at reviewing the fun¬ 
damental mechanisms devised so far, mainly in transition-metal oxides, which 
lie at the origin of multiferroicity. 
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1. Introduction 


Magnetic and ferroelectric materials, characterized by a spontaneous symme¬ 
try breaking that causes the appearance of a switchable long-range magnetic or 
dipolar ordering below a critical temperature, are ubiquitous in modem science 
and present-day technology for their diverse applications, rang ing fr om data stor¬ 


age to magnetic and ferroelectric random-access memories [llil2 


,l3i|40. In recent 


years, the quest for magnetoelectric multifunctional integration within a single 


Email addresses: paolo. baroneSspin. cnr. it (Paolo Barone), 
Silvia. picozziOspin. cnr. it (Silvia Picozzi) 

Preprint submitted to Comptes Rendus - Physique 


February 11, 2015 














material has motivated a renewed interest in the elass of so-ealled multiferroie 


materials, di^laying the simultaneous presenee of two or more spontaneous fer- 
roie phases yD. The intrinsie eombination of magnetism and ferroeleetrieity in 
this elass of eompounds ealls for novel deviee paradigms exploiting their eross- 
eoupled effeet, e.g. allowing to eontrol magnetization (polarization) by an ex¬ 
ternal eleetrie (magnetie) field. The intense researeh aetivity on this elass of 
materials, both from the experimental and theoretieal sides, is testified by the 
number of exeellent review papers on the topie, opening promising prospeets for 
applieations also beyond the field of magnetie and ferroeleetrie materials, e.g. 
for information and energy harvesting teehnologies[6, 7. 8. 9. kI 11. 121. 

From the theoretieal point of view, the intrinsieally eoupled and multifune- 
tional nature of multiferroies poses numerous fundamental ehallenges. In faet, if 
on the one hand the theoretieal understanding of magnetie insulators is rather 
well established, on the other hand a rigorous mieroseopie theory for ferro- 
eleetries has been formulated only in the last few deeades. The mieroseopie 
origin of magnetism is basieally the same in all magnetie insulators, being re¬ 
lated to the presenee of loealized eleetrons, mostly in the partially filled d ox f 
shells of transition-metal or rare-earth ions, whieh have a eorresponding loeal¬ 
ized magnetie moment. Exehange interaetions between the loealized moments, 
usually aeeounted for in general Heisenberg-like spin models, lead to magnetie 
order lllSll . For ferroeleetries, the situation is quite different, and one major difh- 
eulty eomes from the faet that eleetrie dipoles in extended-state systems are well 
defined only in the very speeial ease of a fully ionie material, while in general 
the eleetrie polarization is a global property of the matter that eannot be de- 
eomposed in loealized eontributions. A rigorous definition of ferroeleetrie bulk 
polarization has been provided only in the early 90s via a quantum-meehanieal 


















approach based on Berry phases 1114 
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16 n . whieh may explain the diffieulties 


eneountered in identifying possible mieroseopie meehanisms responsible for fer- 
roeleetrieity. In faet, the mieroseopie origin of ferroeleetrieity is not unique, and 
the eoexistenee and interplay between struetural distorsions and eleetronie de¬ 
grees of freedom (spin, eharge, and orbital) has been proven to play a eentral 
role in the diverse meehanisms devised so far. In this respeet, first-prineiples 
ealeulations based on the density-funetional theory (DFT) have played an im¬ 
portant role in the deseription, understanding as well as predietion of magnetie, 
ferroeleetrie and magnetoeleetrie properties of mutiferroies, due to its ability 
to deseribe the many aetive degrees of freedom within a eomparable level of 
1911 (see also the ehapter of Ghosez and eollaborators). 


aeeuraeyllV, 


IE 


When diseussing the origin and mieroseopie meehanisms for multiferroieity, 
therefore, the main problem lies in the origin of ferroeleetrieity, its eoexistenee 
with and its eoupling to the magnetie ordering. Generally speaking, one ean 
identify two kind of meehanisms leading to the appearenee of ferroeleetrie po¬ 
larization. The first is substantially ionie/displaeive, in the sense that it involves 
lattiee distortions of ions earrying different eharges, where the struetural defor¬ 
mation lies at the basis of the spaee-inversion symmetry breaking. The seeond 
kind of meehanisms, on the other hand, involves primarily eleetronie degrees 
of freedom, whieh are responsible for the spaee-inversion symmetry breaking of 
the eleetronie ground state even in eentrosymmetrie struetures. Sueh a elassifiea- 
tion is to many extents an approximation, sinee the two meehanisms are tightly 
intertwined and a signature of symmetry-lowering in both ionie and eleetronie 
seetors is usually found in the very same material. A somewhat related, but pos¬ 
sibly less arbitrary, distinetion is that between proper and improper ferroeleetrie¬ 
ity: in improper ferroeleetries, the spontaneous polarization ean be eonsidered 



















as a by-product of another structural (or electronic) primary phase transition, as 
opposed to proper ferroelectrics, where the symmetry lowering can be ascribed 
primarily to polar distortions. In the following, we will review the microscopic 
mechanisms for multiferroicity that have been proposed and identified so far, 
paying special attention to the understanding of the origin of ferroelectricity and 
its coexistence/coupling with magnetic ordering. 


2. Microscopic mechanisms for proper multiferroics 


2.1. Conventional ferroelectric perovskites: revisiting the exclusion rule 

Most of the long-known (conventional) ferroelectrics are perovskite transition- 
metal oxides (BaTi03, PbTiOa, KNb03, Pb(ZrTi)03) in which the polar distor¬ 
tion mainly involves an off-centering of the perovskite B-site transition-metal 


cation showing an empty d shell, whic 
valency with surrounding oxvgens OOL 


Cochran [I22I. 


1 is prone to establish some degree of co- 


21II1 . On the basis of a lattice shell-model. 


2311 showed that covalent interactions tend to be short range, while 


the ionic electrostatic interactions are long (in fact infinite) range. The competi¬ 
tion between short-range forces, which tend to favour high-symmetry phases, 
and long-range Coulomb interactions, which destabilize the centrosymmetric 
structure, is then influenced by the strong covalency that softens the cation-0 
repulsion, leading to an offcentering of the transition-metal ion towards one (or 
three) oxygen(s), at the expense of weakening the bonds with other oxygen ions. 
In a simplified local approach, the microscopic origin of the cation offcentering 
can be deduced by estimating the energy associated to the deformation of the 
covalent bonds in a ligand-cation-ligand geometry as 
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( 1 ) 


















where tp^ is the hopping matrix element deseribing hybridization interaetions 
between the eation and the 0-p states, A is the eharge-transfer gap, u deseribes 


the distortion and the ehanges in hybric 


ization are aeeounted for as tpd{u) 


tpd{l ± gu) in the linear approximation 112411 . If the elastie energy eost of the 
lattiee distortion, ~ Ku^/I, is smaller than the quadratie eovaleney energy gain, 
then the loeal distortion is energetieally favourable, whereas the different eharges 
assoeiated to different ions may explain qualitatively the appearenee of loeal 
eleetrie dipoles 3- 

A more refined (but still substantially loeal) approaeh is based on the vi- 
bronie eoupling between ground and exeited eleetronie states, whose expression 
ean be dedueed by general symmetry eonsiderations, that is known to trigger 
(pseudo) Jahn-Teller distortions 11251 12611 . In this approaeh, the adiabatie potential 
energy surfaee near the high-symmetry eonfiguration eontains a yibronie eontri- 
bution K = Kq + Ky. is a positive eontribution that eoherently ineludes the 
loeal odd displaeements of all the atoms, whieh is a long-range (whole erystal) 
feature, whereas the matrix elements of Ky, whieh eontains negligible intereell 
interaetion terms and is found to be always negative, are nonzero when the over¬ 
lap between the wave funetions of the ground state of atoms of one sublattiee 
(oxygens) with the exeited state of the atoms of the other sublattiee (transition- 
metal eations) inereases due to the nuelear displaeements u, thus enhaneing the 
eovaleney. In this framework, the mieroseopie meehanism leading to ferroelee- 


'Due to the covalent character of the bonds, however, a quantitative estimate of bulk P is 
only accessible via the Berry-phase formulation in the framework of the modern theory of po¬ 
larization; in this context, each ion carries dynamical Born effective charges (BECs), instead 
of nominal static ones, and the strong covalent effect is usually reflected in anomalously large 
BECs. 
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trie transition in a crystal is directly related to both its atomic and electronic 
structures, where the origin of polarization appears to be substantially loeal but 
depends strongly on long-range interaetions. 


Besides providing an example of the intrinsieally eoupled nature of ferro- 
eleetrieity in erystal systems, the ease of eonventional perovksites also allows 
to introduee one of the major theoretieal ehallenge that has been faeed in the 
field of magnetoeleetrie multiferroies, i.e. the understanding of whieh eondi- 
tions must be met in order to eombine magnetism and ferroeleetrieity in a sin¬ 
gle homogeneous material. In faet, for a long time the two phenomena have 
been tought to be mutually exelusive and ehemieally ineompatible in perovskite 
oxides, due to the requirement of having empty (partially filled) d states in or¬ 


der to have ferroeleetrieity (magnetism) [127 


280. However, a eareful analysis 


based on the pseudo Jahn-Teller effeet has shown that sueh exelusion rule is 
less striet than expeeted even in the elass of proper ferroeleetrie oxides, where 
some speeifie d” and spin eonfigurations ean allow, on the basis of symmetry- 
allowed (lattiee-eleetron) yibronie eouplings, for the eoexistenee of magnetism 
and ferroeleetrieity [|29L [SOi ISjJ]. A partieularly interesting ease is that of high- 
spin d^ eonfiguration, a situation realized in alkaline-earth manganites sueh as 
CaMn03 or SrMnOa; these systems were shown to possess a weak ferroeleetrie 
instability mediated by a eovaleney-driven meehanism, whieh at ambient eondi- 
tions is hindered by other energetieally favourable distorsions, namely nonpolar 
rotation/tiltings of BOg oetahedra |132i[SSi [3^ . Applying strain or ehemieal pres¬ 
sure eould in prineiple tune the balanee between different struetural instabilities. 


allowing for ferroeleetrieity to emerge, as inc 


CaMnOs films [l35h and bulk Sri_vBa^MnO^ [l36l] . Sinee both magnetie and ferro 


eed reeently shown for strained 


eleetrie instabilities are related to the same manganese B eations, the magneto- 
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electric coupling is expected to be strong in these svstems OTL 13811 . even though 
antiferromagnetic phases are likely to appear; on the other hand, critical temper¬ 
atures may be far larger than in frustrated magnets (e.g. a Tc ~ 185 ^ has been 
reported for Sro.sBao.sMnOaDS^). 


2.2. Ferroelectricity due to lone pairs 

The recent boost in the multiferroic field can be reasonably ascribed also 
to the successful realization of BiFeO^ thin films with enhanced multiferroic 


properties by the Ramesh’s group ll39ll . This perovskite oxide, one of the most 


celebrated multiferroic materials due both to its very large ferroelectric polar¬ 
ization and to its high ferroelectric as well as magnetic critical temperatures, 
apparently escaped from the aforementioned exclusion rule. However, it turned 
out that the microscopic origin of ferroelectricity in this system is completely 
different, not involving polar offcentering of transition-metal B-site cations; in¬ 
stead, the main instability was related to the presence of stereochemically active 


lone pairs on the A-site Bi ions 1140 


411]. The ferroelectricity arises here from on¬ 


site sp rehybridization of the two 6s electrons of bismuth that do not participate 
in chemical bonds, thus showing a high polarizability whose particular orien¬ 
tation may create local dipoles and trigger the onset of the observed very large 
ferroelectric polarization (P ~ 100 pCIcm HA'M ). It is interesting to note that 
the presence of nonpolar tilting/rotation of the BOe octahedra, that ultimately 
seems to prevent the appearence of ferroelectricity in orthorhombic alkali-metal 
manganites, is not competing with the ferroelectric instability in BiFeOs, where 
polar distortions involve primarly A-ions; on the contrary, a primary role of the 
BOe rotational patterns has been recently suggested in the ferroelectric phase of 
this compound ll43ll . On the other hand, magnetism is guaranteed by the B-site 
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magnetic (d^) ions. Due to the independent origin of ferroelectricity and 
magnetism in this type of multiferroics, one could have expected the magneto¬ 
electric coupling to be rather small; however, a number of recent observations 
indicate that the magnetoelectric coupling in BiFeOa can be significant and lead 
to unusual and very interesting effects, suggesting also an important role played 


by the large ferroelectric 


of the magnetic ordering [144 


£0 


arization in determining a long-period modulation 




m 


£ZD. 


Beside BiFe 03 , the lone-pair mechanism could be realized in other related 


systems. Indeed, it was first theoretically proposed for BiMnOafldSD, even though 


later computational and experimental structural optimizations revealed that the 
ferroelectric phase is not stable and suggested that BiMnOa belongs to the cen- 


trosymmetric C2/c structured. 


son . The ferroelectric (or antiferroelectric) na¬ 


ture of this material is still debated, as several details - such as oxygen stochiom- 
etry - seem to pla y a crucial role in determining both the structure and the prop¬ 


erties of BiMnO^ llS ill . Recently, it has been also proposed that BiMnOa could 


behave as an improper ferroelectric, where ferroelectricity could develop from a 
specific antiferromagnetic ordering which would break the inversion symmetry 
even in a centrosymmetric crystal structure 152. On the other hand, a similar 


mechanism 


has 


PbNiOaiSd, 


55, 


Deen suggested to be at play in the relatively new multiferroic 


56, 


5711 . despite the fact that the nominal valence of Pb is 4+, 


thus no lone pairs should be active on A-ions; DFT-based results have however 
clarified that the Pb formal valence is perturbed by Pb 65 - O 2p hybridization 
which results in enhanced ferroelectric instability and predicted large polariza¬ 
tion (P ~ 100/iC/cm^). 
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3. Microscopic mechanisms for improper multiferroics 


3.1. Hybrid improper ferroelectricity 

In the quest for general meehanisms whieh may give rise to both ferroelee- 
trie and magnetie orderings indueed by the same lattiee instability, a signifieant 
theoretieal advanee eame from a eareful reeonsideration of rotational distortions 
in transition-metal ABO 3 perovskites. Oetahedron rotations are in faet ubiq¬ 
uitous in perovskites and related materials, and are known to strongly couple 
to magnetic properties, because they buckle the inter-octahedral B-O-B bond 
angles which mediate the interplay of electronic (spin and orbital) degrees of 


freedom 115 811 ■ Although they are usually found to compete with ferroelectric in¬ 


stabilities in conventional ferroelectric oxides, it has been recently shown that 
some combination of such rotations can induce local polar displacements of A- 
site cations through effective force fields exerted by O ions, even though the 
distortional patterns that are most frequently found in simple perovskites result 


in cooperative antipolar (antiferrodistortive) distortions 11591] . However, this ob¬ 


servation suggests the exploration of related systems where the effect could be 
exploited to engineer an improper ferroelectric state arising from suitable com¬ 
binations of proper rotational distortions (which are individually nonpolar). The 
key ingredient is a trilinear coupling between a polar mode P (that individually 
would be energetically unfavourable) with two unstable nonpolar lattice modes 
Ri^Ri with different symmetries, whose combined symmetry properties Ri ©i ?2 
are the same of the polar one. The mechanism, that was originally proposed to 


explain ferroelectricity in a Aurivillius layered compound llbOl] . has been called 


“hybrid improper ferroelectricity”, highlighting the improper origin of ferroelec¬ 
tric polarization and the role of coupled nonpolar distortions 116in . The fact that 

the polarization is proportional to the product of two distortions implies that it 
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will flip its direction if either or R 2 , though not both, is reversed. Interest¬ 
ingly, the trilinear coupling may also serve as the interlink between magnetic 
and ferroelectric properties, i.e. provide a mechanism for linear magnetoelectric 
effects. This has been thoroughly discussed in Ref. 1^ . where a layered per- 
ovskite Ca 3 Mn 207 served as a model system to illustrate the coupling between 
ferroelectricity. The theoretical magnetic ground state displays a canted anti¬ 
ferromagnetic configuration with a weak ferromagnetic moment, arising from a 
relativistic Dzyaloshinskii-Moriya interaction that is strongly dependent on the 
buckling of Mn-O-Mn angles. Interestingly, a linear magnetoelectric effect was 
predicted to appear, mediated by the rotational modes which are responsible at 
the same time for the ferroelectric polarization and the (weak) ferromagnetic 
moment. Since these nonpolar distortions usually manifest near or above room 
temperature, there appears to be no fundamental limitation on the temperature 
range over which this mechanism could give rise to multiferroic order, thus mo¬ 
tivating an intense research activity aiming at identifying suitable design rules to 
single out promising candidates, as discussed in the chapter by Ghosez et al. 

The mechanism has been identified in several materials that are potential im 


proper ferroelectric and multiferroic, including perovskite artificial superlattices 1162. 


layered (AB 03 ) 2(;^)1 


A/A’ cation orderings 1164 


6111 or double AA’BB’Og perovskites with specific 


6511 . Interestingly, it has been recently shown how a 


very similar mechanism can be active at twin walls naturally occurring in many 


ABO 3 centre symmetric perovksites (such as CaTi 03 or PbZrQ^f ll66ll . where po¬ 


lar displacements have been predicted [l67ll and then observed to occur at the 


wallrtbS 


6911 . Due to the structural origin of such effect, a wall polarization is 


expected to develop also at twin walls of magnetic nonpolar oxides, such as 

CaMn 03 or LaFe 03 , that would be possibly coupled to the magnetic ordering. 
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Figure 1: Hybrid improper ferroelectricity arising from a cooperative interplay of rotation and 
tilting distortions as usually found in the Pnma orthorhombic structure, (a) The distorted oxygen 
positions are responsible for a net force field acting on A-site ions (asterisks mark the oxygens 
that get closest to the central A-cation). In the Pnma structure, characterized by an pattern 

in Glazer’s notation, the A ions undergo antipolar displacements, which usually do not give rise 
to a net electric polarization. When such antipolar displacements are not compensated, as in the 
case of ordered layers of A and A' ions (b) or in the layered structure found in the Ruddlesden- 
Popper phases (c), a bulk polarization can develop from the interplay of the nonpolar rotational 
distortions. 








































Eventually, other nonpolar lattiee modes eould be ineluded beside rotations and 
tilting; as a matter of faet, hybrid improper ferroeleetrieity and magnetoeleetrie 
effeets - in prineiple allowing for an eleetrieal eontrol of magnetization - have 
been predieted in related metal-organie frameworks with perovskite strueture, 
where the nonpolar mode mediating both effeets is a Jahn-Teller distortion of 
oxygen oetahedra eoupled with a tilting mode|^ 71]. As opposed to their in- 
organie eounterparts, metal-organie frameworks offer the possibility of varying 
their additional degrees of freedom, arising from organie/inorganie duality, and 
are eurrently the objeet of an inereasing interest as potential novel multiferroie 
materials ll72l] . 


3.2. Magnetically induced ferroeleetrieity 

Among the elass of improper multiferroies, materials displaying magnetieally 
indueed ferroeleetrie polarization have been objeet to an intense researeh aetiv- 


ity in the last deeade, due to their large intrinsie magnetoeleetrie eouplingC, 
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IJ, ITj]. Generally speaking, the ferroeleetrie polarization in these materi¬ 


als arises as an improper (or seeondary) transition indueed by a speeifie mag- 
netie ordering, whieh is itself responsible for the symmetry lowering and loss 
of the spaee inversion even on top of a eentrosymmetrie strueture. The way by 
whieh the magnetie ordering breaks the inversion symmerty is two-fold, sinee 
the symmetry breaking may oeeur in the magnetie eleetronie ground state (even 
in eentrosymmetrie erystal lattiees) or may be eaused by struetural distortions 
indueed by speeifie eonfigurations of the magnetie moments (a meehanism gen¬ 
erally ealled exehange strietion). 

At the basis of the first kind of meehanisms lies the intrinsie eharge anisotropy 
of individual d orbital states, that is ultimately responsible for both eovalent ef- 
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fects and exchange interactions in transition-metal oxides, whereas the coupling 
between spin and orbital degrees of freedom can be mediated by the atomic spin- 
orbit coupling or by correlation effects such as the Hund coupling. In this respect, 
model approaches have been devised to identify possible mechanisms by which 
a local electric dipole may develop from electronic spin degrees of freedom, typ¬ 
ically in a cluster model taking into account the transition-metal orbital states 
and their hybridization with surrounding ligand ions, mediating the effective d-d 


hopping processes and crystal-field splittings 1174 




77 


78 


m, 


80 , 


MD- Un¬ 


der general assumptions (substantially related to the time-reversal invariance), it 
is possible to show that the electric dipole of a spin dimer with inversion sym¬ 
metry can be always written in a cluster expansion approach as: 


P - Pl{Si) ■+■ Pr{Sr) + PlrSl X Sr, 


( 2 ) 


where the local contributions must be a quadratic form of the corresponding 


local spin, P^, 


aaSa/s, a = I, T and Of,/? labeling the cartesian com¬ 


ponents x. V. 7 117811 ■ while /?/,• is a 3 x 3 matrix. Assuming some microscopic 


model for the cluster, the explicit dependence of the coefficients jia and pir 
on microscopic quantities can be derived. The easiest way is to consider two 
transition-metal ions interacting via some ligand, described by a model Hamil¬ 
tonian H = Hl Ht Hsoc + Htm + Hu, where Hu (Htm) describes the lo¬ 
cal p (d) states on the ligand (transition-metal) sites, H, accounts for the elec¬ 
tron transfer through the ligands, which is described by an hybridization ma¬ 
trix depending on the p, d orbitals involved and on their relative position, while 
Hsoc describes the spin-orbit coupling on the transition-metal ions. The last 
term Hu = -U J^ri=ir^a • describes a local Zeeman-like term which allows 


to tune the spin-dimer configuration. In a three-ion linear cluster with the lig- 
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and occupying the middle of the bond, two eontributions to polarization have 


been predieted 174, 


75L 


76l. I 77 II . arising from the mixing of d orbitals indueed 


by the spin-orbit eoupling and hybridization elfeets mediated by the ligand oxy¬ 
gen. This mixing triggers the onset of a loeal dipole in the eluster, intimately 
eonneeted to the shapes of involved d and p orbitals. From the nonzero overlap 
integral {da,yz\y\Pz), a transverse eomponent of the loeal eleetrie dipole is found, 


P oc e X S/ X Sr, 


(3) 


where e is the unit veetor parallel to the bond direetion l74l l76h. Following Kat- 
sura, Nagaosa and Balatsky, the eleetrie dipole is said to be indueed by a spin- 
eurrent meehanism, sinee the veetor produet S/xS^ is proportional to the spin eur- 
rent ji, where the Dzyaloshinskii-Moriya veetor d aets as its veetor potential ly . 
Within the same eluster geometry, a purely eleetronie longitudinal eomponent of 
the eleetrie dipole also arises from the nonzero overlap integrals {da^xy\^\py) and 
{da,zx\^\Pz)- This eontribution to P results from the variation of the d-p hybridiza¬ 
tion mediated by the spin-orbit eoupling, with a predieted funetional form in the 


linear three-ion olusterflTSi 


m 


P oc (S/ ■ e) Si - {Sr ■ e) Sr 


(4) 


showing its loeal eharaeter (involving only single-spin terms) as opposed to the 
previously diseussed meehanisms for inter-site spin-dependent hybridization ef- 
feet. From general symmetry argument, it ean be shown that the only non-zero 
eoeffieients entering in Eq. ^ are = (Ci,0,0), nf’ = 7r|^ = (0,Ci/2,0), 
ttJ' = jfj^ = (0,0, Cl/2) and pir{2,3) = -pir{3,2) = C 2 in the linear three-ion 


model, and Eq. [2]reduees to the former equations (|3]), (141) 117 811 . However, the 

explieit expressions of the eoeffieients pir and iia depend on the speeifie geom- 
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etry of the cluster and/or the material under considerations, affecting local elec¬ 
tronic structure, crystal-field splittings and hopping interactions. For instance, 
according to the conventional spin-current mechanism of Eq. ([3]), the helical 


spin-spiral configuration in triangular-lattice antiferromagnetic CuFeOo[l82l. 


m 


or MnE llSdIl should not support any ferroelectric polarization. On the other hand, 
a combined DFT and cluster model approach has shown that indeed a general¬ 
ized spin-current mechanism can be active in the specific triangular-lattice ge¬ 
ometry. In fact, the presence of two ligands connecting the magnetic ions is 
responsible for the appearence of other nonzero terms in the pir matrix, allowing 


for the appearence of a polarization even in the helical spin configuration llVSh . 


On the other hand, the local contribution to P arising from the spin-dependent 
hybridization mechanism in tetrahedral crystal fields, as found in melilite crys¬ 
tals, has been shown to behave as P = Jja(S ■ Caf Ca, where a = 1, ..,4 labels 
four oxygens surrounding the central magnetic ion and Ca are unit vectors par¬ 


allel to the metal-oxygen bonds [l 851 


8011 . whereas a more complex dependence 


have been found in Cu 20 Se 03 , displaying CuOj trigonal bipyramids and square 


pyramids [l 8 6h . 


Even though a unified cluster-based approach that attempts to describe both 


pure electronic and ion-displacement contributions has been recently proposed [l79n . 


magnetostrictive effects have been usually described in the general context of 
spin-model Hamiltonians, where the effect of structural distortions on exchange 
coupling constants has to be suitably included, depending on the material un¬ 
der consideration. A necessary condition for magnetically-induced polar dis¬ 
tortions is that the magnetic ordering itself breaks the inversion symmetry, i.e. 
only inhomogenous or frustrated spin systems can support a polar displacement. 
Microscopically, lattice distortions in the presence of competing magnetic inter- 
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Figure 2: (a) and (b) Two possible choices of cluster models describing a spin dimer with one or 
two bridging ligand ions. If inversion symmetry is preserved, two contributions to local electric 
dipole are expected, Pi depending only on the local spin configuration on each transition-metal 
site and Pir generally being proportional to Si x Sr- In a linear three-ion model (a), the intra¬ 
site contribution gives rise to a purely longitudinal contribution to polarization, whereas the spin 
current induces a transverse local dipole. In the four-site cluster (b), the longitudinal P, devel¬ 
ops in the plane of the cluster, its direction depending on the spin-dimer configuration; on the 
other hand, the spin-current polarization is not always perpendicular to the spin current, having 
different in-plane and out-of-plane components, (c) Magnetostrictive effect in the T T i i Ising 
model, (d), (e) Magnetostrictive effects in Mn-O chains representative of R MnOa, with strong 
oxygen buckling distortions. The spins on each Mn ion experience a superexchange coupling 
J ^ Jq + > 0, which is invariant under inversion symmetry and depends only on even 

powers of z„; assuming z„ - (-l)"zo + small O displacements are energetically favourable if 
6zn - COS On, where 6„ is the angle between neighbouring spins and | 6z^ is the 

elastic energy associated to the distortion. When considering the antisymmetric Dzyaloshinskii- 
Moriya interaction, the staggering of the Dzyaloshinskii vector along the chain implies that the 
corresponding magnetic energy is minimized by 6zn {AJk) sin0„; in this case, a collective O 
shift is allowed only by a spiral spin configuration, since the vector product x 5„+i has the 
same sign for all pairs of neighboring spins 
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actions can then appear in order to maximize the gain in magnetie energy. A 
simple example is provided by an Ising-spin chain with competing ferromag- 
netie nearest-neighbour and antiferromagnetie next-nearest-neighbour interae- 
tions, where an inhomogeneous T T i i spin eonfiguration ean be realized. If 
some energy is gained, e.g., by the shortening of a ferromagnetie bond, the ehain 
will distort in sueh a way that bonds between parallel spins shorten and those 
between antiparallel ones streteh, causing polar displacements that would result 
in a net bulk polarization if different eharges are carried by eaeh site along the 
ehain. The simplest realization of this magnetostrictive meehanism can be found 
in CaaCoMnOe, where Co^^^ and magnetie eations alternate along one¬ 


dimensional ehainslISV. 


8811 . On the other hand, the (super)exehange proeesses 


in transition-metal oxides are mediated by the ligand oxygens bridging the metal 
eations, and they generally depend on ionie relative positions. Sinee oxygens 
earry nominal negative charge (as opposed to the positive one earried by mag¬ 
netic ions) their colleetive shift can induce a finite polarization. Different mag¬ 
netic configurations can be responsible for different polar displacements through 
the symmetrie or antisymmetrie exehange interaetion, as diseussed heuristieally 
in Ref. 1891] for rare-earth manganites RMnOa. In order to maximize the mag¬ 
netie energy gain, the Mn-O-Mn bond angle mediating the exehange interaetions 
may be inereased or deereased, eausing a polar displaeements of oxygens. A 
possible spin ordering leading to polar displacement is the t T i i eonfiguration, 
where the Mn-O-Mn angle increases (deereases) at ferro- or antiferromagnetie 
alternating bonds, leading to a colleetive displacement of ligand O driven by 
the symmetric exchange interaetion. On the other hand, the antisymmetrie (or 
Dzyaloshinskii-Moriya) exehange interaetion is responsible for a polar magne- 

tostrietion if a eyeloidal spin eonfiguration sets in in the orthorhombie strueture 
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of rare-earth manganites, as shown in Fig. [2l The polarization driven by such 
inverse Dzyaloshinskii-Moriya mechanism appears to have the same dependence 
on spin configuration as the spin-current polarization derived before; however, 
the latter has been derived in a frozen ionic configuration and describes the purely 
elecronic contribution to the ferroelectric polarization, whereas the former de¬ 
scribes a magnetostrictive mechanism where ionic displacements are induced 
mainly by the magnetic energy optimization. In principle, both an electronic and 
an ionic mechanism can contribute to the total polarization, even though it is not 
easy to predict which contribution is larger. 


3.3. Ferroelectricity due to charge/orbital ordering 

Improper ferroelectricity is also in principle allowed in materials displaying 
specific charge or orbital orderings which break the inversion symmetry, indi¬ 
vidually or in combination with a magnetic ordering. Charge ordering is a rather 
ubiquitous effect in transition-metal oxides, being often observed in materials 
with ions having formally a mixed valence like Fe or Ni. It is also likely to appear 
in low-dimensional charge-transfer organic salts. Generally speaking, different 
forms of charge ordering may appear. Typically in materials with rather local¬ 
ized electrons, the electronic charge disproportionation takes place mainly on top 
of ions, leading to a so-called site-centered charge ordering. On the other hand, 
the Peierls transition often observed, e.g., in quasi-one-dimensional organic salts 
is characterized by a bond-centered charge-density wave, implying that charge 
disproportionation takes place roughly on alternating strong and weak bonds be¬ 
tween ions. Each time a given charge-ordered configuration breaks the space- 
inversion symmetry — for instance, in the presence of a mixed site-centered and 
bond-centered charge-density wave, as proposed in [|%] — a macroscopic po- 
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larization is likely to appear 19111 . There are, however, no general rules for the 
realization of either one situation or the other (or a eombination of the two), and 
eharge ordering effeets have to be eonsidered separately ease by ease. The sim¬ 
plest situation is realized when eharge-ordered ions oeeupy sites belonging to a 
strueture that laeks inversion symmetry, as in the ease of some of the proposed 
low-temperature struetures of magnetite Fe^QA [l92ll or in the double-layer strue¬ 


ture of LuFe 204 , displaying Fe^'^-rieh and Fe^'^-rieh layers II93L 19411 . When tak¬ 
ing into aeeount the eovaleney, sueh ionie pieture ean be modified and a purely 
eleetronie eontribution to polarization ean appear; this effeet ean be understood 
in the framework of the Berry-phase theory of polarization by looking at the 
Wannier-funetion eenters, deseribing the eenter of eharge of the eontinuous elee¬ 
tronie eharge density, whieh in the presenee of eovalent effeets may move from 


eentrosymmetrie positions 


appearenee of loeal dipoles [l95h . 


site-eentered or bond-eentered), thus leading to the 


Similarly, orbital ordering may also eause the appearenee of eleetronie eon¬ 
tribution to polarization when the eleetronie eharge density aequires a noneen- 
trosymmetrie distribution that ean be also traeed by the positions of the eor- 
responding Wannier-funetion eenters 11951] . In its simplest realization, orbital- 
indueed polarization ean appear when an orbital ordering develops in a system 
with bond dimerization and inequivalent hopping integrals; the orbital-ordered 
state ean then be viewed as a superposition of two inequivalent eharge-loealization 
phenomena in different orbital seetors, as shown in Fig. [3 ]d), leading to uneom- 


pensated loeal eleetronie dipoles whieh sum up to a bulk polarization 11981] . This 


meehanism eould be realized, e.g., in vanadate spinels sueh as CdV 204 , showing 
both Peierls-like bond dimerizations and an ordering of dy^, d^x orbitals. A dif 


ferent seenario has been proposed for undoped and half-doped manganites[l96 

19 


















Figure 3; (a) and (b) Sketch of a neutral one-dimensional chain displaying a mixed site-bond 
charge-ordering wave (a) or an orbital ordering in the presence of bond-dimerization (b). In both 
cases the Wannier-function centers are displaced with respect to inversion-symmetric positions, 
highlighted by vertical dotted lines, thus giving rise to uncompensated local electric dipoles and, 
hence, to a macroscopic polarization, (c) Schematic representation of electronic contribution to 
polarization expected in the E-type antiferromagnetic configuration found in some undoped man- 
ganites, such as HoMn03; electronic hopping is constrained by strong Hund’s rule only along 
one-dimensional zig-zag chains in MnOa layers, and is further asymmetrized by the onset of 
orbital ordering triggered by Jahn-Teller and correlation effects. The Cg-like Wannier-function 
centers are then coherently displaced from the Mn sites, giving rise to a Berry-phase polariza¬ 


tion parallel to the chain direction ll96l] . (d) Schematic diagram of the electronic polarization 


and buckling-induced mechanism predicted in the CE-type antiferromagnetic configuration ob¬ 
served in half-doped manganites, such as Pi'o.sCao.sMnOa; the interplay of correlation effects and 
buckling distortion is responsible for orbital ordering and a bond-dimerization along each ferro¬ 
magnetic chain, where the Mn-O-Mn angle between Mn ions with parallel occupied orbitals is 
increased more than the angle between Mn ions with perpendicular orbitals. Polar displacements 
of bridging oxygen cause the appearance of an ionic P perpendicular to the chain, while orbital 
ordering leads to an electronic Berry-phase P parallel to the chain direction ll97l] . 
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9711 . whose electronic properties can be properly described in the limit of 


infinite Hund coupling, which induces an infinite intra-atomic splitting between 
minority- and majority-spin states. The relevant low-energy model is a degener¬ 
ate double-exchange Hamiltonian describing eg spinless electrons whose motion 
is determined by the underlying magnetic configuration of the almost localized 
t 2 g spins. The combination of E-type antiferromagnetic ordering and Jahn-Teller 
lattice distortions, as found in undoped manganites, remove the orbital degen¬ 
eracy leading to an ordered state with alternating d^x'^-r^ /along zig-zag 
chains of parallel spins in each Mn02 layer. In turn, this leads to a change of 
the Berry phase associated to the alternating right-handed and left-handed mo¬ 
tion of the Cg electrons around each site displaying the orbital orderin g, w hich 


is ultimately related to the appearence of an electronic polarization [l9a 


990. As 


for the case of spin-induced electronic polarization, the mechanism relies on the 
coupling between spin and orbital degrees of freedom, where the magnetic con¬ 
figuration determines the way in which d states hybridize and electrons hop in 
the lattice and a local dipole develops depending on the shape of individual or¬ 
bitals and their hybridization through oxygen’s p orbitals; in this case, however, 
the strength of such coupling is determined by Hund’s interaction, rather than 
the weaker spin-orbit coupling, with a consequent larger value of t he expecte d 


contribution to polarization (as indeed predicted by DFT calculations lIlOOl 


imD). 


4. Conclusions 

In this chapter, we have presented a tutorial introduction to the microscopic 
mechanisms of multiferroicity that have been unveiled so far. The recent theoret¬ 
ical advances in the field seem to suggest that the multiferroic phenomenology 

is much more common than what originally believed. The reason probably lies 
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in the variety of mechanisms leading to ferroelectricity and to its coexistence 
with other kind of orderings, involving also electronic degrees of freedom. Even 
though several theoretical approaches have been adopted to identify the origin of 
multiferroicity in its diverse forms, the strength of the microscopic mechanisms 
and their efficiency in leading to large polarization can be most likely assessed 
by resorting to first-principles techniques, as thoroughly discussed in the chapter 
by Ghosez and collaborators. Nonetheless, the intuitive understanding hope¬ 
fully provided in this chapter might provide the basic knowledge required for the 
search and optimization of multiferroic properties in novel materials beyond the 
class of perovskite oxides, including even more complex oxides, such as mate¬ 
rials with spinel or melilite structure, or organic-inorganic hybrid systems as the 
recently proposed metal-organic frameworks. 
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